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The super symmetric intertwining relations with second order supercharges 
allow to investigate new two-dimensional model which is not amenable to 
standard separation of variables. The corresponding potential being the 
two-dimensional generalization of well known one-dimensional Poschl- Teller 
model is proven to be exactly solvable for arbitrary integer value of parameter 
p : all its bound state energy eigenvalues are found analytically, and the 
algorithm for analytical calculation of all wave functions is given. The shape 
invariance of the model and its integrability are of essential importance to 
Qh' obtain these results. 
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■ 1 Introduction. 



The beautiful idea of supersymmetry (SUSY) was first introduced pQ and developed in 
Quantum Field Theory and Elementary Particle Theory at the seventies of the last century. 
During these years supersymmetry became one of the most popular and promising branches 
of modern High Energy Physics [2]. 

Supersymmetry was also studied in the simplest toy model of (0 + 1) Quantum Field The- 
ory (i.e. in nonrelativistic Quantum Mechanics) in order to clarify some delicate problems 
^ ■ of spontaneous supersymmetry breaking [3]. Very soon, this by-product of supersymmet- 
rical Quantum Field Theory became a new independent tool to study many problems in 
Quantum Mechanics itself [3J. In particular, the notions of SUSY intertwining relations [I], 
[5] and of shape invariance [6] provided both new methods to derive some old results and 
to obtain new interesting results. As an example, all previously known one-dimensional 
exactly solvable potentials were reproduced as potentials obeying the shape invariance [7]. 
In its turn, SUSY intertwining relations were successfully used [8], [H], [ID], [EE] in two- 
dimensional Quantum Mechanics to obtain a variety of partially ( quasi-exact ly) solvable^ 
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models which are not amenable to conventional separation of variables [13]. Right up to re- 
cent time, the latter method was the sole practical tool to solve analytically two-dimensional 
(and higher-dimensional) quantum problems. 

Thus, supersymmetrical approach can be considered as a new method [9], \±0\ to solve (at 
least, partially) some two-dimensional Schrodinger equations. The procedure can be called as 
SUSY separation of variables. It is realizable for the models where equation for zero modes of 
second order supercharges allows separation of variables [10]. More of that, after separation 
of variables one-dimensional equations must be exactly solvable too. Another procedure of 
SUSY separation of variables works if one of the partner Hamiltonians does allow standard 
separation of variables due to special choice of parameters. Then, SUSY intertwining rela- 
tions may allow to obtain eigenfunctions of the second partner Hamiltonian, which is not 
amenable to standard separation of variables. This approach was used successfully [H] for 
the two-dimensional generalization of Morse potential with integer or half-integer values of 
parameter. 

The present paper provides new exactly solvable two-dimensional model with potential 
depending on three parameters, one of which has to be integer. To solve the problem, 
it will be necessary to explore essentially both main ingredients of SUSY Quantum Me- 
chanics: SUSY intertwining relations and shape invariance. Schematically, to solve the 
Schrodinger equation with potential V(x; A, B,p) depending on parameters A,B,p three 
steps will be done. First, to find such exclusive value of parameter (actually, p = 1,) that 
initial Hamiltonian H(p = 1) does allow conventional separation of variables. Second, using 
SUSY intertwining relations and shape invariance, to build eigenfunctions for Hamiltonians 
H(p), p = 2,3, which are not already amenable to separation of variables. And finally, 
to prove that all constructed wave functions are normalizable and that no extra levels exist. 
The structure of the paper is the following. In Section 2, the model and its main properties 
are formulated, and the scheme of investigation is reviewed. The separation of variables 
for first Hamiltonian H(p = 1) is performed and some delicate properties of potential are 
discussed in Section 3. The zero modes of the supercharge are built in Section 4, and they 
are used for construction of wave functions in Section 5. Normalizability of wave functions 
is studied in Section 6, where the absence of any other bound states was proven. A few 
examples of wave functions for low values of parameter p are given in Section 7. Conclusions 
includes the comparison of obtained results with the limiting case which is explicitly solvable. 
Rather cumbersome calculation of coefficients necessary for wave functions and spectrum are 
presented in Appendix. 

and corresponding wave functions are known. Such models take up an intermediate place between exactly 
solvable ones and models with unknown spectra |12j . 
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2 Formulation of the Model and the General Scheme. 

We consider the intertwining relations of the form 

Q-H = HQ-- HQ + = Q + H, (1) 

where H and H are two-dimensional Hamiltonians of the Schrodinger type 

H = -(d\ + d 2 2 ) + V(x u x 2 )- H = -(dl + d 2 2 ) + V(x u x 2 ), (2) 

and intertwining operators Q ± are second-order differential operators. A number of models 
of this kind were investigated in the series of papers [8], [9], [10], [JTJ. In [H], it was proven 
that one of them - the generalized two-dimensional Morse - possesses exact solvability. Here 
we shall show that one more model [JT] involved in intertwining relations (CQ) is also exactly 
solvable. It reads: 

H(p) = -dl - dl- 2p(p - 1) (cosh~ 2 (a; + ) + cosh" 2 (a;_)) + 

ki (sinh _2 (x 2 ) — cosh~ 2 (xi)) + k 2 (cosh _2 (x 2 ) — sinh~ 2 (a; 1 )) , (3) 

H(p) = -dl - dl- 2p(p + 1) (cosh" 2 (x+) + cosh _2 (a;_)) + 

k\ (sinh~ 2 (a;2) — cosh~ 2 (xi)) + k 2 (cosh~ 2 (x 2 ) — sinh _2 (a;i)) , (4) 

Q ± = dl -dl±2p (tanh(x+) + tanh(a;_)) d^± 

2p (tanh(x_) — tanh(x + )) d 2 + 4p 2 tanh(x + ) tanh(x_) + 

k\ (sinh _2 (x2) + cosh~ 2 (xi)) + k 2 (cosh _2 (x2) + sinh~ 2 (xi)) , (5) 

where p and k± j2 are real parameters, so far arbitrary, and x± = X\ ± x 2 . The Hamiltonians 
in ©, © can be represented in the form: 

H P _ T (xi) + H P _ T (x 2 ) + f(x 1 ,x 2 ), (6) 

where Hp_t(x) are well known one-dimensional Poschl- Teller Hamiltonians, and f(xi,x 2 ) 
- specific term mixing X\ and x 2 variables in potentials. Due to this expansion, potentials 
V(xi,x 2 ), V(xi,x 2 ) may be considered [TT] as a two-dimensional generalization of Poschl- 
Teller potential. These models are shape-invariant [B] with respect to the parameter p: 

H(p+l) = H(p). (7) 

Properties ([TJ) and ([7]) will be essential for the proof of exact solvability of the model for 
positive integer values of p. 
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We remark that, by construction, these models are integrable since from intertwining 
relations ([1]) it follows that 

[H,R}=0, [H,R] = 0; R = Q + Q~, R = Q~Q + , (8) 

with symmetry operators of fourth order in momenta. 

The general scheme to determine the spectrum of the model ([3]) could be adopted from 
the paper |14| . where the full spectrum of two-dimensional generalization of Morse poten- 
tial [9], [TO] was found. In the present context, the plan of construction could be the fol- 
lowing: we start with H(p = 1) and find all normalizable solutions ^(x;p = 1) for the 
corresponding Schrodinger equation (Section 3) as far as it is amenable to separation of vari- 
ables. Then, by means of intertwining relations ([I]), we find eigenfunctions ^(x;p = 1) 
of H(l). In general, they might be of two types [§]: some of them are inherited from 
H(l) as: ^(x; 1) = Q~(l)^f(x; 1), and others are zero modes of the intertwining opera- 
tor: Q + (l)^(x; 1) = 0. In such a way we obtain all eigenfunctions of H{2). Due to the 
shape-invariance ([7]) of the model, ^(x;p + 1) = ^/(x;p), and therefore, we have calcu- 
lated already 2). Following this strategy step by step, we expect to find the eigenfunc- 
tions and eigenvalues for the Hamiltonians H(p) = H(p + 1) with arbitrary integer values 
p = 1, 2, .... At each step, the full variety of eigenfunctions of H(p + 1) will belong to one of 
two classes: 1) each normalizable wave function ty(x;p) leads to normalizable wave function 
ty(x;p+ 1) = ^(x;p) = Q~ (p)^/(x;p); 2) the same Hamiltonian H(p+ 1) has also some 
number of extra normalizable functions which are specific linear combinations of zero modes 
Q(x;p) of the operator Q + {p). We shall see below that this plan has to be modified suitably 
for the case of our present model, but the main ideas will be analogous to that of [9J, [10J, 

D~Q- 



3 Separation of Variables for H(p = 1). 

For the Hamiltonian H(l) the standard procedure of separation of variables in Carte- 
sian coordinates can be applied. Looking for the solutions of the Schrodinger equation 
H(l)fy(x; 1) = E^(x; 1) in the form 1) = f](xi)p(x2), one obtains two one-dimensional 
equations for unknown functions p, 77 

-7/'(xi)-( — — + — -J — )vM = £77(2:1); (9) 
\cosh X\ sinn x\ J 

-p"(x 2 )+( — -| — + — — ) p(x 2 ) = sp(x 2 ), (10) 
\cosh £2 sinn x 2 J 

where prime denotes the derivative of the function with respect to its argument, e + e = E 
is the energy value for H(l), and both e and e must be negative for the discrete part of 
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the spectrum. Thus, we have to consider solutions of one- dimensional Scrodinger equations 
(|9|)- (|T0|) with Poschl- Teller potentials Vp_r(x). It is convenient to replace parameters ki, k 2 
by A, B according to k\ = B(B — 1); k 2 = — A(A — 1). Avoiding the case of fall onto center 
[T5] . we shall restrict ourselves with reasonably attracting singularity in Vp-t(xi), Vp-t(x 2 ) 
with coefficients ki G (-1/4,0), k 2 G (0,1/4), i.e. it is sufficient to take A,B G (0,1/2). 
The substitution r](xi) = sinh A (xi) cosh' B (xi)F(xi), and the subsequent change of variable 
Xi to z = — sinh 2 (a;i), turns ([9]) into the hypergeometric equation for the function F(z) : 



z(l 



z ffl + (x + i_^ + B + 1 ) z )^ + (_i ( x + ^_i £)FW = o. 



dz 2 



The pair of independent solutions for the given value of e reads (see 2.3.1(1) in |T6]): 



T]f\x) = sinh A (x) cosh i( (a;) 2 -Fi 



A + 5 + A + B 



— ]A+ -; -sinh 2 (x; 



'ir 



?7^ 2 ' ) (x) = sinh 1 A (a;) cosh 5 (a;) 2-P1 



1 - A + 5 



£ 1- A + B 



-e 3 



— A; — sinh (x] 
(12) 

The similarity of expressions (JTTJ and ( fT2l) reflects the obvious symmetry of potential in 
OH]) under A — >■ (1 — A). The potential under consideration obeys also the similar symmetry 
under — > (1 — B). But the corresponding independent solutions are related to solutions 
fTTTj) . (|12p according to relations between hypergeometric functions (see 2.1.4(23) in [T6]). 
The formulae analogous to (iTTj) and (fl2|) hold also for p- 1 ^^), but with the necessary 
changes A ^ B, B A and e — >• e. 

To provide the normalizability of \I/(x; 1), both p and 77 must be normalizable. To ana- 
lyze the possible bound states of H{jp — 1) it will be sufficient to consider the asymptotic 
behaviour for large \x%\. 

The general solution r] £ (x) is a linear combination: 



r) e (x) = axrjP(x) + a 2 r)¥ ) (x) 



(2), 



(13) 



with arbitrary constants a±,a 2 . The asymptotic behaviour of the analytic continuation of 
hypergeometric functions for large z = — sinh 2 a; (see 2.10(2), 2.10(5) in [16]) reads: 



2 F 1 (a, b; c; z) = B x (a, b, c) 



-z)- a +0{{-z 



,-a-l\ 



-B 2 (a, b, c) 



{-z)- b +0{{-z)-^) 



(14) 



where constants B\ 2 are expressed in terms of Gamma functions: 

T(c)T(b-a) T(c)T(a-b) 
Bi(a, b, c) = r] B 2 {a, b, c) 



T(b)T(c-a) 



T(a)T(c~b) 
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Substitution of (JHj) into (fTTj) . (fl2|) and (|T3l) gives asymptotically two groups of terms in 
r] £ (x) : proportional to (— z)~^~^l 2 {l + 0(;z -1 )) and proportional to (— z) + ^~^^ 2 (l + 0(z^ 1 )) 1 
correspondingly. In order to forbid the growing term in wave function, one has to require 
the coefficient to vanish: 

D (A + B + yf=i A + B- yf=i . 1 
° A { 2 ' 2 ' + 2 

+« 2 5 2 ( - 2 , - 2 , - 2 ~A)=0. 

In general, there are two options to fulfil this requirement: 

n D fA + B + ^E A + B- yf=i , 1\ n 

^ = 0; B 2 ( , ! 2- A J =0; 

In the case of arbitrary A, B, these conditions can be achieved by means of suitable choices 
of energy values e, due to Gamma functions in denominators of coefficient B 2 : arguments a 
or (c — b) of these Gamma functions must be equal — n with n = 0, 1.... Just this condition 
might give the energies of bound states. But one can check easily, that in our present case of 
parameters A, B 6 (0, 1/2), these conditions can not be fulfilled for positive values of \f^e. 

Thus, we are not able to kill the growing terms in asymptotic of rj e and p§. Therefore, the 
Hamiltonian H(p = 1) has no bound states because of asymptotic behaviour at large \x\, and 
the first source for construction of eigenfunctions of H(l) : ty(x; 1) = Q~(l)^/(x; 1) - does 
not work. We stress that this statement depends crucially on a chosen region for values of 
A, B, which in its turn was dictated by conditions on both potentials Vx, V 2 , simultaneously. 
Taken separately, these Hamiltonians would have bound states, but for different values of 
A,B. 

We notice also that the conclusion above does not depend on the behaviour of solutions 
at the singular point x = 0. Nevertheless, we will discuss the x — > asymptotic of rfP'^ 1 
here, since it will be necessary for the analysis in subsequent Sections. The point is that 
(in contrast to standard situation of nonsingular potentials) both solutions (llljl . ( I12j) have 
zero limit at the origin x — >■ for A e (0,1/2). Namely, their behaviour is r/^ ~ x A 
and T]^ ~ x l ~ A . This is the typical situation of the so called "limit circle" kind (see [T7] . 
Appendix to Section 10.1), which was widely discussed in the literature in the context of 
one-dimensional potential (the so-called Calogero potential) g/x 2 on the semi-axis or on 
the whole axis (e.g., see [TS], [IS], [20] )• In such a case, there is continuous freedom in 
choosing (among many opportunities) some "good kind of behaviour" for wave functions. 
The resulting spectrum of the model depends on this choice, thereby defining the kind of 
its quantization. The preferences are usually motivated by physical arguments [19], |18j . 



6 



[2Tj . In the case of Hamiltonian H(p = 1) this problem is of little importance due to 
asymptotic behaviour of solutions at infinity discussed above. One more remark concerns 
the extension of solutions to negative semiaxis: it is reasonable to choose an odd way, taking 
r](— \x\) = — rj(\x\). This choice provides continuity of the derivative rj'(x) at the origin. 

We have to remark that from mathematical point of view, both one-dimensional 
Schrodinger operator with Poschl- Teller potential in fl9]), f JTOj) and two-dimensional oper- 
ator in ([3]), (j4|) produce rather nontrivial problem. It is possible to check that both of them 
(in two-dimensional case, due to Green's identity of vector calculus), are symmetric opera- 
tors, but in strictly mathematical approach, they are unbounded and not self-adjoint for the 
conventional choice of smooth functions with a compact support (dense in L2) as a domain 
D(H P ^ T (x)). Similarly to the analysis given in [19] for the case V = a/x 2 , the self-adjoint 
extension of Hp_j<(x) includes also the functions from L 2 with specific asymptotic at the 
singular point x = 0. For details, we refer readers to the papers [19], [18], [20] and references 
therein, where one will find also the description of some paradoxes induced by too naive 
approach to singular potentials of a/x 2 type. 

4 Construction of Zero Modes of Q + . 

As it was mentioned above, the second possible source of eigenfunctions for H(p) are the 
zero modes of operator Q + (p). It is well known [9], [10], that the subspace of zero modes of 
Q + (p) is invariant under the action of H(p). This means that if Qk{x;p) is the zero-mode of 
Q + (p), i.e. Q + (p)flk(x; p) = 0, then due to intertwining relations (JTJ) 

TV 

H(p)Q k (x;p) = y^C fc A(f;p) (15) 

(Cki is the (N + 1) x (N + 1) matrix with complex elements). If the matrix can be 
diagonalized by some matrix B : 

BC = AB; A = diag(E , E 1 , E N ), (16) 

the functions 

N 

$i(x; p) = J2 B iMx; p) (17) 

k=0 

are the eigenfunctions of H(p) : H(p)^i(x;p) = Ei^i(x;p). 

At first, one needs to calculate Qi(x;p). For this purpose, it is useful to perform the 
similarity transformation, which will help to separate variables: 

q +( p ) = e Px(i)Q+( p y-px(i). Q i = e- px ^Ui, (18) 
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q + (p)ui{x;p) 







(19) 



where q + reads 

q + = d\ - d\ + k x (sinh~ 2 (x 2 ) + cosrT 2 (a;i)) + k 2 (cosh~ 2 (x 2 ) + sinh" 2 (xi)) . (20) 

The choice of the function \ provides that (fl9l) is amenable to separation of variables 
in Cartesian coordinates. The two-dimensional equation ( Fl9l is equivalent to the pair of 
one-dimensional ones if one takes Ui = r](xi)p(x2), and it appears that they are exactly the 
equations fl9|)- f|T0|) . but with e = e. The solutions can be written as linear combinations of 



where r\ £ (and analogously, p £ ) must be built from solutions (jTTj) . (fI2|) . Of course, we are 
interested only in normalizable zero modes Q(x;p) of the two-dimensional operator Q + (p), 
but the normalizability condition, in comparison with Section 3, is essentially less restrictive 
now: 



The factor exp (— 2px{x)) in (122]) is exponentially decreasing at infinity in all directions 
on the plane, and it is able to compensate even growing functions u £ . Due to asymptotic 
equivalence coshx ~ sinhx at infinity, asymptotical behaviour of the integrand of ( 1221) can 
be represented as 



Therefore, the functions Q are normalizable for arbitrary values of p and e, satisfying: e > 
—p 2 . This fact has to be taken into account in calculation of the spectrum of H(p + 1) (see 
Section 5). 

But at first, we must define the variety of functions Q £ , which may be used for construction 
of actual wave functions. In this context, functions f] £ {xi) and p e (x 2 ) are the auxiliary objects 
for construction of zero modes Q according to fflBl) . (T2TT) . Therefore, all four possible combi- 
nations can be used, in general. The first of them ile(x) = exp ( — px(x) \ri £ l \xi)p ( i\x2) 



Us = r]e(,xi)p £ (x 2 ), 



(21) 




(22) 
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is: 



rif^(x) == ± (cosh(x + ) cosh(x_)) p sinh A (|xi|) cosh B (|a;i| 



■^mh B {\x 2 \)cosh A {\x 2 \) 2 F 1 [a £ ,b E -A+^,z 1 j - 2 F 1 L £ ,b £ ;B + ^;z 2 ) (23) 



where 



and the sign ± depends on a quarter on a plane (x 1 ,x 2 ), according to the choice at the end 
of Section 3. Other zero modes Q £ 2 \x), Q £ 3 \x), Q^\x) are obtained from fljp(x) by means 

of substitutions of pairs of parameters ( 1 — A, B ] for Qs(x), of ( A, 1 — B ] for Q,^\x), 



and ^1 - A, 1 - 5 J for ^ 4) (x), instead of (^L, Bj in ^(f). 

5 Eigenfunctions of the Hamiltonian H. 

In this Section we shall look for linear combinations of zero modes of Q + {p), which are simul- 
taneously the eigenfunctions of the Hamiltonian H(p) = H(p + 1) in 01]). Being interested in 
the discrete energy spectrum E n , we suppose that the corresponding wave functions are built 
from the finite number of zero modes £l £k \x;p); 7 = 1,2,3,4, with parameters at = a £k in 
f )24|) . numbering by discrete values k = 0,1, 2,..., so that the constants are ordered 
as: ao> a\> ... > a N (~,). We suppose also that four kinds of such wave functions exist: each 
is built from the corresponding zero modes with fixed value of 7 (the value of 7 defines 
behaviour at the origin). According to (IT5|) . 

N 

H(p)^(x) = J2 C ^W> ( 25 ) 

where c£j are constants, and N also depend on 7 = 1, 2, 3, 4. 

Performing with H(p) the similarity transformation analogous to that with Q + (p) in (j!8p . 
one obtains: 

h(p) = 4*MH(p)e-*xM = -dl -dl + D- -^ - -J* + _?* + -^ - Ap\ 

cosh x\ sinli x\ cosh x 2 sinn x 2 

where the mixing operator D is defined as: 

2v 

D = — (smh(2zi )dx + sinh(2a; 2 )<9 2 ). (26) 

cosh cosh x_ 
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Then, exploring (|9]), ( fTUl) . the action of h(p) on (p) (see its definition in ( 1181) and (121 

can be expressed as: 

h(p)uW(p) = 2^2p(A + B)+e k - 2p 2 ^£W + sinh A (xi) cosh B (xi) sinh B (x 2 ) • 
• cosh A (x 2 )i)^ 2 Fi(a fc , b k ; A + 1/2; - sinh 2 x x ) 2 Fi(a k , b k ; B + 1/2; - sinh 2 x 2 ) 
and (T25|) takes the form: 

2^2p(A + B) + e fc - 2p 2 ^) 2 F l (a k , b k ; A + 1/2; z x ) 2 F 1 (a k , b k ; B + 1/2; z 2 ) + 



+D ^ 2 Fi(a k , b k ; A + 1/2; z x ) 2 F x {a kl b k ; B + 1/2; z 2 ) 

N 

= Yl C ki 2 i? i(«i> ^ + 1/2; *i) 2^1(0*, 6t; 5 + 1/2; z 2 ). 

i=0 

After straightforward calculations, ([25]) can be rewritten as: 

2^2p(A + B) +e k - 2p 2 ^ 2 Fx{a k ,b k ; c x , z x ) ■ 2 F x (a k , b k ; c 2 ; z 2 ) + 
+T^k (^r 1 2Fi(a k + l,b k + 1; t* + 1; zi) • 2 Fi(a fc> 6 fc ; c 2 ; z 2 ) + 
2 Fi(a fc + 1, b k + 1; c 2 + 1; z 2 ) • 2 ^i(afc, &fc; ci; Zi) 

^A r ^(1) 



Z 2 (l-^2) 
C2 



Z)i=o G fci 2^iK,6i;ci;^i) • 2 i 5 i(oi, 6;; c 2 ; z 2 ), (27) 



where 



* = a + \,* = b + \, a, B ^ = d±g+^ , (a 



For z 2 = 027]) reads: 

2 ^2p(A + 5) - 2p 2 + ^ aFiCafc, 6 fe ; c i; Zl ) + *9&B Zl ■ 2 Fi(a fc + 1, b k + 1; ci + 1; z x ) = 

= ^Cj»-*F 1 (a i ,b i Wz 1 ). (29) 
In the z\ — oo limit, the largest power in the l.h.s. of ( 1291) is (— z\)~ ak . Therefore, 

CkJ = for i > k, (30) 
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i.e. Eq.( l29|) takes the form: 

2^2p(A + B) -2p 2 + e k ^j ■ 2 F 1 (a k , b k ; a; z x ) + ^^ Zl ■ 2 F x (a k + 1, b k + 1; c x + 1; z x ) = 

= ELo C 1I ) ■ 2F 1 (a i ,b i ;c 1 ]z 1 ). (31) 
Further, for particular values z x = and k = 0, it gives: 

Cfi } = 2 (2p(4 + £) - 2p 2 + so) . (32) 



Substitution of ( 132]) back into (l3Til with arbitrary £i and = leads to ao • &o = 0. Since 
all 6 fe are positive, the only opportunity is a = 0. Comparing next powers in Eq.(!3Tl) for 
Zx — > — oo, we obtain: 

a k + 1 = a k -i- 

Together with a = 0, this relation uniquely defines all values of a k : 

a k = -k; k = 0,l,...,iV (1) . (33) 

In turn, comparison of coefficients of (— zi)~ ak in Eq. (l3"Ti) . gives values of elements C^ kk . 
Due to fl30|) . matrix is triangular. Its diagonal elements cj^jj coincide with elements 
of diagonal matrix A in (fIB"j) . and therefore, C^j} gives a part of the eigenvalues of discrete 
energy spectrum of the Hamiltonians H (p) = H(p + 1) : 

E k \v) = E i k 1 \p+l) = -2{{A + B + 2k-pf +p 2 y, fc = 0,l,...,iV (1) . (34) 

It is clear from (1341) that the lowest energy state corresponds to the maximal k, i.e. to 
k = which can be defined from conditions of normalizability of Q^J(x). These conditions 
were formulated in Section 4, and they can be rewritten now as: 



fcW < \ p -A-B)- AT« 
2 



\ip-A-B) 



(35) 



where [c] means the integer part of c. 

Analogously one can construct three other kinds of energy levels E ( k l] (p) of H(p) = 
H(p + 1) by replacing everywhere above (A,B) by (1 — A, B), or by (A, 1 — B), or by 
(1 — A, 1 — B), correspondingly. The result is the following: 

Ef\p) = E ( k 2 \p+l) = -2{[l-A + B + 2k-pf+p 2 y, fc( 2 > = 0,l,...,jV< 2 >; (36) 

Efip) = E^\p + 1) = -2^(1 + A- B + 2k-p) 2 +p 2 y, k^ = 0,l,..-,N^ ; (37) 

E { k A \p) = E { k A \p+l) = -2^2-A-B + 2k-p) 2 +p 2 y, feW = 0, 1, JV^, (38) 
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where: 

A; (2) < ^(p-l + A-B); (39) 
k {3) < -(p-l-A + B); (40) 

k^ < ^(p-2 + A + B). (41) 

The energy spectra and the corresponding wave functions for several lowest values of p will 
be given in Section 7. 

According to (jTTl) . the eigenfunctions of H(p + 1): 

N 

^\x;p + 1) = ¥^(x;p) = Y,B^\x-p)- 7 = 1,2,3,4 (42) 

can be written explicitly only after calculation of coefficients C|j in (1251) and, after that, of 
-B^ from ffTBT) . The first step - calculation of C^? with general k, i - is given in Appendix, 
while the calculation of seems to be rather complicated in a general form. Instead, this 
will be done explicitly for small values of p in Section 7. 

Thus, each Hamiltonian H(p + 1) (for > 0) possesses + + + + 4 
bound states ^\x;p + 1) with energy levels E^\p+ 1), k = 0, 1, iV (7) , 7 = 1,2,3,4, 
which are absent in the spectrum of H(p). As we know, due to SUSY intertwining relations 
([T]) (see also Section 2), each of these wave functions produce the tower of extra eigenfunctions 
for higher Hamiltonians H(p + n + 1), n = 1,2, ... with the same energy values E% (p +1). 
These wave functions ^^(x;p + n + l) are built by the action of n operators Q~ (p + m), m = 
l,2,..,n: 

¥^(x;p + n + l) = ¥^(x;p + n) = Q-(p + n)Q-(p + n-l)...Q-(p+l)^\x;p+l), (43) 

and their indices indicate the number k among bound states of original Hamiltonian 
H(p + 1), and the number n of one after another acting operators Q~ . 

6 Normalizability of the wave functions. 

According to results of previous Sections, the Hamiltonian H(p + 1) = H(p) has two classes 
of bound state wave functions. The second one (in terminology of Section 2) ^^\x;p + 
1), k = 0, 1, is produced by normalizable zero modes of Q + via their suitable linear 
combinations. The first class is obtained from the eigenfunctions of lower Hamiltonians by 
means of operators Q~ . In notations introduced above, they are: 

= Q~(P)Q~(P-1)---Q~{n + I)*tt&n + 1), (44) 
n = l,2,...,p-l, 
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with energy values Em (n + 1) (see (1341) ). The restrictions for values of m depend on 7 : 

m < (n - A - B)/2, 7 = 1; m < (n - 1 + A - B)/2, 7 = 2; 
m < (n - 1 - A + B)/2, 7 = 3; m < (n - 2 + A + B), 7 = 4. 

This is an appropriate point to remark that the situation of general position corresponds to 
the simple spectrum of H(p + 1), which consists of levels E^\p + 1) (see fl34|) - fl38|) ) for 
(p+ 1) and levels Em\n + 1) for the states (T4"4"|) . But nothing prohibits from the possible 
occasional degeneracy of the spectrum for some specific values of parameters. Indeed, this 
situation is nongeneric: an occasional degeneracy of levels may occur only for some single 
values of parameters A, B. In such a case, the degeneracy can be removed easily by an 
arbitrary small variations of A, B. 

The normalizability of functions of the second class is obvious by construction, but this 
property for the wave functions (jUJ) will be proven now. The Hamiltonian H(p) has the 
symmetry operator R(p) = Q~ (p)Q + {p) (see ([8])), and in turn, H(p+ 1)— its own symmetry 
operator R{p + 1) = Q + (p + l)Q~(p + 1). As far as these Hamiltonians coincide (shape 
invariance) H(p+ 1) = H(p), the corresponding symmetry operators must coincide as well, 
but up to the function of the Hamiltonian. Indeed, by straightforward calculation one obtains 
the relation: 



R(p + 1) - R{p) = 8(2p + 1) ( H{p) + 2{2p 2 + 2p + 1) 



(45) 



which will help to analyze the normalizability. 

The norm of the arbitrary wave function V&^l _ n \(x;p + 1) (T4"4"|) can be written as: 



•Q-(p)g-(p-l)...Q-(7i + 2)g-(n + l)*W(f;n + l)). 



To simplify it, one may explore Eq.( l45|) and its consequence: 




Q + (n + l)Q"(n + 1) - Q-{n)Q + {n) )Q-(n)Q~(n - 1) . . . Q"(m) 



Q-{n)Q-{n-l)...Q-{m)Y 



where V 



run 



is the function of the Hamiltonian: 



r 



run 
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The following relation can be derived by induction: 




(if; p + 1) of second class 



for H(p + 1) are proportional to the norms of wave functions \l/m (if; u + 1) of the first class 
for the Hamiltonian H(n + 1) : 



It is easy to check explicitly that the coefficient of proportionality is positive. Hence, as 
far as the initial state ^fm (n + 1) is normalizable by the construction, any wave function 
(p-n) (P + -*-) * s n o rm alizable too. 

One more statement is necessary to prove in order to be sure that the full variety of 
eigenfunctions for H(p + 1) was constructed above. Namely, we must prove that no addi- 
tional normalizable wave functions exist besides those in f|4*2|) . f l4"3j) . Starting from the lowest 
Hamiltonians, let us suppose that H{2) has such additional eigenfunction $(2), which differs 
from the linear combination of zero modes of Q + (l). Then, it follows from the intertwin- 
ing relations, that Q + (l)Q(2) must satisfy the Schrodinger equation with Hamiltonian H(l). 
One may check that supercharges Q (p) do not change the normalizability neither at infinity, 
nor at coordinate axes Xi,x% : the detail analysis is presented below in the next paragraphs 
of this Section. As we already know from Section 3, the Hamiltonian H(l) has no bound 
states at all, and therefore, our supposition was wrong. Let us suppose now that the first 
Hamiltonian possessing such additional state $(p + 1) is H(p+ 1), while all previous Hamil- 
tonians H(p), H{p — 1), ...H(2) have bound states of the forms (H2J), (]4"3l) . only. Then, due 
to intertwining relations, Q + {p)Q{p + 1) is the eigenfunction ^(p) of H(p), and therefore, 
coincides either with ^k{p) or with ty^i(p), by our assumption. For simplicity, we do not 
consider here the case of possible degeneracy of levels of H(p) (the conclusion will be the 
same in this case). Acting by Q + (p) onto $(p + 1) and using the relation (|4"5l) . one obtains 
by straightforward calculations that for both options, ^(p) is proportional to Q + ^ n ,i{p + 1) 
with some suitable n. Therefore, the wave function <3>(p + 1) coincides with ^> n i(p + 1) up 
to zero modes of Q + {p) : 




(x; p + 1) || 2 = 6A(2n + 1) ( E$ (n + 1) + 2n 2 + 2n + 1 




q=n+l 



$(f;p + 1) = ci(p + l)*„,i (z;p + 1) + c 2 (p + l)* fc (f;p + 1), 
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where c\^{p + 1) are constants. 

Thus, the problem is reduced to the question: whether the operators Q ± {p) are able 
to change the normalizability of functions. If they are not, no additional normalizable 
eigenfunctions of H(p + 1) exist. It is evident from the explicit expressions (jSJ) of Q (p) and 
from taking into account the exponential decreasing of \l/ at infinity, that Q ± can not violate 
integrability of |\I/| 2 at ±00. 

More difficult problem arises in the neighborhood of x% — y and/or x 2 — * 0. The part of 
Q ± linear in derivatives coincides with the operator D, defined in ff26l) . In the limit x 2 — > 0, 
x\ ^ 0, it is: 

D ~ 4p(tanhxi)(9i — 4p(l — tanh 2 Xi)x2d2-, 

i.e. it does not change the asymptotic behaviour of the function. Analogous conclusion is 
true in the limit 22 — ► 0, X\ 7^ 0. To analyze the limit when both x\^ 0, it is convenient 
to use the polar coordinates x + = Rcosip,x_ = Rsinip. Asymptotically, D for R — » is: 



D ~ 4p(wa.(2<p)Rd R + cos(2v?)<9, 




i.e. in this limit D can not change the behaviour of function as well. 

Coming back to the operators Q ± , the only terms which could in principle change the 
behaviour of function at xi — > and/or x 2 ~ > are: 



Q ± (p) ~ - (^-dl - k 2 sinh- 2 (xi)^ + (-d$ + h sinh- 2 (. 
ese terms with (1H1). (HI), we observe the same narts fait 



Comparing these terms with ([3]), (jl]), we observe the same parts (although with different 
signs) in asymptotical expressions: 

H(p) ~ (^-9 2 - fc 2 sinh- 2 (x 1 )^ + (^-d 2 2 + fc lS inh- 2 (a;2 

Since ty(x;p) are eigenfunctions of H(p), Q ± (p) are not able to change the behaviour of 
and the absence of any additional wave functions besides that of (1421) . ( 143|) types was thus 
proven. 



7 Examples. 

The explicit expressions for matrix elements from (|16p . which are necessary to build the 
eigenfunctions f|42l) of seem to be difficult to present in a general form. Nevertheless, 

the problem can be solved straightforwardly for low values of p. By means of separation of 
variables, we demonstrated in Section 3, that the Hamiltonian H{1) has no bound states. 

The next Hamiltonian H (2) (it corresponds to p = 1 in formulas above) has two bound 
states: one bound state with A;^ = 0, due to inequality (!35j) with p = 1, and the second 
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bound state with k^ = or A;( 3 ) = 0, due to inequalities (|39|) . ( HOjl . depending on the 
positivity of (A — B) or (B — A). Of course, these bound states are of the second class, i.e. 
are built from the zero modes: 

¥o\x; 2) ~ Q ( q \x; 1) = ± (cosh(x+) cosh(a;_)) _1 ■ 

■ sinh A (|a;i|) cosh B (|a;i|) sinh B (|x 2 |) cosh A (|x 2 |) (46) 

with energy E^(2) = -2 ((A + B - l) 2 + l^j , and (for A > B) 

*j, 2) (x; 2) ~ fi[, 2) (f; 1) = ± (cosh(x+) cosh(a;_)) _1 ■ 

• sinh^dxil) cosh 5 ( | x x |) sinh B (|x 2 |) cosh^flafcl) (47) 



with energy Ef } {2) = -2\[B - A) 2 + lj . No bound states of the first class exist in this 
case. 

For the next value p = 2, i.e. for the Hamiltonian H(3), four bound states are of the 
second class being built by the zero modes 

¥< 7) (f;3)~n< 7) (f;2); 7 = 1,2,3,4 

with energies 

(3) = -2 ( {A + B - 2f + lY £<; 2) (3) = -2 f (5 - A - l) 2 + 1 



E® (3) = -2^4-5-l) 2 + lJ; eJ 4) (3) = -2 ^(A + 5) 2 + 1 

But in this case, two wave functions of the first class also can be built from (146]) and ( 14"?]) by 
the procedure (j4"3l : 

^ 1 ) (f;3)~g-(2)vI/W(x;2); 7 = 1,2 A > B. (48) 

Their energies E^l\3) coincide with E^\2); 7 = 1,2 above. 

The Hamiltonian H(p + 1) with p = 3 for A > B has six states which are built 
from bound states of H(3) by means of operator Q~(4). They have the same energies 
E^(3), £j 2) (3), E®(3), ^J 4) (3) and E^{2), E^ ] {2). As for the second class bound states, six 
such bound states exist: k^ = 0, 1; k^ = 0, 1; k^ = k^ = 0. This set includes two wave 
functions coinciding with Qq(x; 4); Q^\x; 4), and four other wave functions have to be built 
as linear combinations of pairs of zero modes ^ 1} (x;4), ^(x;^ and ft[, 2) (x;4), ^S 2) (f;4), 
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since iV (1) = iV (2) = 1 for p 
defined by (J55J), §S]): 



4. The matrix elements of 2 x 2 triangular matrix are 



(i) 

00 



E«(4) 



(A+B-3) 2 +9 



Cg } = ^(4) = -2 



,(1) 



(A+£-l) 2 +9 



C 



(i) 

10 



-24. 



These elements are necessary to determine coefficients 5^ for (|42p . From Eq. (1161) . one can 
find that for cff ^ Cff, as in the present case, the matrix is also triangular B§± = 0. 
More of that, while 



B 



(i) 

10 



a 



(i) 

10 



(1) 



a 



— B y 

l) ii 



(l) 
oo 



the last coefficient i?^ is arbitrary. This fact is not discouraging, since the linear combination 



2J) for ^q 1 '' (4) includes only one term (with arbitrary B [ 



by unity norm of ^q 1 ^). The second linear combination (l4"2l) for ty\ l> (4:) includes both B 



(i> 
oo 



Thus, the value of Bqq is fixed 



and , which are proportional to each other. The absolute values of these coefficients will 
also be fixed by normalization of the wave function. Analogous calculations can be easily 
repeated for 7 = 2. 
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8 Conclusions. 

An exhaustive procedure of analytical solution of two-dimensional generalization of Poschl- 
Teller model with integer values of parameter p was presented above. Being based on SUSY 
intertwining relations and shape invariance of the model, the procedure replaces the standard 
method of separation of variables which is not applicable here, and it can be considered as 
a special - SUSY - separation of variables. 

In order to confirm obtained results for H(p + 1), it is useful to compare them with the 
limiting case which possesses the direct solution by means of separation of variables. Indeed, 
if the parameters A, B, which originally belong to the interval (0, 1/2), are chosen on the 
limit of range A, B — y 0, the procedure above (starting from Sect. 3) does not work. But 
due to conventional separation of variables, the Hamiltonian H(p + 1) from (J3j) is reduced 
(up to a trivial multiplier 2) to a sum of two one- dimensional Hamiltonians with well known 
reflectionless potentials in variables x± : 

h(p+ l)(x) = -a 2 — • 

cosh x 

The spectra of these Hamiltonians are well known: for p 6 [L, L + 1) they have exactly L 
bound states. To compare the properties of spectra with that in Section 7, one has to explore 
the original inequalities ( 135|) . ( 139|) - ( |4T|) and ( 134|) . ( 136|) - ( 138|) where A and B are written 
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explicitly. The point is that some of bound states described in Section 7 disappear in the 
limit A,B -> 0. 

Thus, H (p + 1) with p — 1 has one bound state with energy E^(2) = -4 in the limiting 
case, since only (1351) can be satisfied with k^ = 0. Taking into account the multiplier 2 
mentioned above, this value of energy coincides with double value of eigenvalue of h(2). 

For p = 2, three bound states of the second class for H(p + 1) exist in the limiting case, 
with jfeW = = £(») = and energies e£\3) = -16; E^ 2) (3) = E^ 3) {3) = -10, and one 
bound state of the first class with energy E^ 1] (2) = -4. The one-dimensional Hamiltonian 
h(3) has two bound states, leading just to four possible bilinear combinations in the standard 
separation of variables. 

For p = 3 with A > B the limiting Hamiltonian H(4) has nine bound states: five of 
second class k^ = 0, 1; k^ = k^ = k^ = 0, and additionally four bound states of first 
class inherited from four bound states of H(3). As it should be, this number coincides with 
3x3 = 9 possible combinations of one-dimensional wave functions. 

Finally, it is necessary to note that the results of the paper might have several destina- 
tions. First, to realize new pure analytical methods of analysis of two-dimensional Quantum 
Mechanics, which are few in number to present day. Second, to use these new methods for 
quantum design in different applications, such as quantum dots, modern nanodevices and 
some cosmological models. Third, successful use of one-dimensional Poschl- Teller poten- 
tial for description of interaction in diatomic molecules signals about perspectives to use its 
two-dimensional analogues in quantum chemistry. Fourth, the complete solvability of the 
present model gives the opportunity to check the validity of different approximate schemes 
in many-particle quantum physics. 
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Appendix: Calculation of the Coefficients Cm- 

We will calculate the coefficients , but with 7 = 2,3,4 can be easily obtained 
by replacing in formulas below A — y (1 — A) etc. The calculation will be started from 
Eq. (l2"Tl) . where parameters are defined according to (12HL (1331 . ( 1281) . It is useful to express 



18 



the hypergeometric functions (with a k = —k) in terms of Jacobi polynomials: 

k\ 

(«Ti) fc 



2 F 1 (a k} b k ;c l ;z l ) = - — 'yrrr-Pk(yi), 



2 Fi(a fe + 1, b k + 1; Cl + 1; z x ) = i*^p£*W>( yi ), 

{a + z) k 

2 Fi(a k ,b k ;c 2 ;z 2 ) = . ' P^' a) (y 2 ), 
[P + J-Jfc 

2 F x (a fe + 1 A + 1; c 2 + 1; * 2 ) = g^pV* 1 *** (y 2 ) , 

\p + Z) k 

where Pokhgammer symbols (r)*. are defined as (T) k = T(T + 1) ■ ...(r + k — 1), and new 
variables and new parameters will be more suitable: 



01,2=1-2*1,2, a = A- p = B--, 



and now: 



b k = 1 + a + (3 + k, ci = 1 + a, c 2 = l + /3, ^e k = 1 + a + (3 + 2k. 
Taking into account, that (a + l)(a + 2)^-1 = (a + 1)^, the l.h.s. of ( 1271) becomes: 

M = ? - (2^(1 + a + (3)-(l + a + P + 2k) 2 - 2p 2 )p^ ) (y 1 ) P ( k M (y 2 ) 



(a + l) fc 03 + l) fc 

- 4p(1 + " + /3 + fc) (d - « w) (k0 ^foO + (i - y 2 2 ) « M+1) (s/ 2 ) Pr\m)) )-(49) 

l/i + 1/2 V /J 

Let us use the relation for Jacobi polynomials 22.17.15 from (22] with n replaced by n + 1 
and /3 by /3 — 1. Multiplying it by (1 + x) and using the relation 22.7.16 from [22], one obtains: 

q + /3 + l\ 2\ r>( a ~t~ii/3"4~i) /™\ 2(?t-~)~c3^(t?-~)~^3^ . 



2 J 2n + a + p 

2w(w + l) (a>/9) , , / jn(n + a) n(n + 1 + g) \ (a>j8) 
2( n + i) + a + /3 ^+ 1 [X) + \ 2 n + a + P 2(n + l) + a + (3j n lj ' 1 j 

Let us write the same relation but for x ==>■ y and a •<=>■ (3, and add it to the initial Eq. (j50p . 
Then, rewriting P n+ i as a combination of P n -\ and P n (according to 22.7.1 from [22J), we 
obtain: 

(1 - J)F^\x)Pt?*\ V ) + (1 - f)P^° +l \y)P^\x) = 

n + l + a + pi 2n + a + (3 \ 
-n(x + y)P^(x)P^ a \y).y 
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Due to this relation, (SSI) takes the form: 



M = , , , r , { -2((l + a + (3 + 2k- P r + p 2 ) P^\yi)Pt a \y 2 ) - 



{a + l) k (P + l) h 
16p(k + a)(fc + 0) 



^ (^foO^fo) + i&W^foi)) }• (51) 



2A; + a + (3 y 1 + 
The r.h.s. of ( JBTT) includes the combination: 

which (by means of recurrent formula 22.7.1 from [22J) satisfy: 

= (w + yi) a -^Pt-f(y.)P^\yi) - ^^-^ (52) 

Ol(fc-l) Ol(fc-l) 

where the following definitions were introduced: 

a ln = 2(n + l)(n + a + /3 + l)(2n + a + 

a 2n = (2n + a + P + l)(a 2 - /3 2 ); 

a 3n = (2n + a + /3)(2n + a + (3 + l)(2n + a + (3 + 2); 

a 4n = 2(n + a)(n + /3)(2n + a + /3 + 2). 



Since 



$o,i = (yi + y 2 )— P (a ^(yi)^' a) (y2), (53) 

aio 



Eqs.(IS2D, §3J give: 

Ol(fe-2) Ol(fc-l) 

+ a 3(fc - 3) a 4(fc _ 1} m + ,,. ( _ 1)fc -i^ fl ^-i)"-^ P (^) (yi)Po ^) (yii) Y(54) 

ai(fc-3) ai(fc-i) «i(fc-2) «io ai(fc_i)...an / 

The coefficients 

,fc + i_ia3i a 4(fc-l) a 4(fc-2)---04(i+l) , _ a 3(fe-l) 



"J- J , 0fc,fc-l = 



Oli Oi(fc-l)dl(A,.-2)--- ( 3'l(i+l) Ol(fc-l) 

allow to write (151]) more compactly: 

fc-i 



i=0 
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Finally, substituting it into ( 149]) . we obtain from ( )27|) the general expressions for desired 
coefficients : 

Cfi? = -2 ((1 + a + /? + 2k - p) 2 + p 2 ) , (55) 

r (i) 16(fc!) 2 p(fc + a)(fc + g)(l + a),(l + 

fc ' i<fc " (*!) 2 (l + a) fc (l + /3) fc (2A; + a + /3) ^ ^ 
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